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ABSTRACT 

We report on the March- April 1997 BeppoSAX observations of Aql X-l, the first to monitor the evolution of the 
spectral and time variability properties of a neutron star soft X-ray transient from the outburst decay to quiescence. 
We observed a fast X-ray flux decay, which brought the source luminosity from ~ 10 36 to ~ 10 33 ergs" 1 in less 
than 10 days. The X-ray spectrum showed a power law high energy tail with photon index T ~ 2 which hardened 
to r ~ 1 - 1.5 as the source reached quiescence. These observations, together with the detection by RossiXTE of 
a periodicity of a few milliseconds during an X-ray burst, likely indicate that the rapid flux decay is caused by 
the onset of the propeller effect arising from the very fast rotation of the neutron star magnetosphere. The X-ray 
luminosity and hard spectrum that characterise the quiescent emission can be consistently interpreted as shock 
emission by a turned-on rotation-powered pulsar. 

Subject headings: stars: neutron — stars: individual (Aql X-l) — pulsars: general — X-ray: stars 



1. INTRODUCTION 

Soft X-Ray Transients (SXRTs), when in outburst, show 
properties similar to those of persistent Low Mass X-Ray Bi- 
naries containing a neutron star (LMXRBs; White et al. 1984; 
Tanaka & Shibazaki 1996; Campana et al. 1998). The large 
variations in the accretion rate that are characteristic of SXRTs 
allow the investigation of a variety of regimes for the neu- 
tron stars in these systems which are inaccessible to persistent 
LMXRBs. While it is clear that, when in outbursts, SXRTs 
are powered by accretion, the origin of the low luminosity X- 
ray emission that has been detected in the quiescent state of 
several SXRTs is still unclear. An interesting possibility is 
that a millisecond radio pulsar (MSP) turns on in the quies- 
cent state of SXRTs (Stella et al. 1994). This would pro- 
vide a "missing link" between persistent LMXRBs and recy- 
cled MSPs. Aql X-l is the most active SXRT known: more 
than 30 X-ray and/or optical outbursts have been detected so 
far. The companion star has been identified with the Kl V vari- 
able star VI 333 Aql and an orbital period of 19 hr has been 
measured (Chevalier and Ilovaisky 1991). The outbursts of Aql 
X-l are generally characterised by a fast rise (5-10 d) followed 
by a slow decay, with an e-folding time of 30-70 d (see Tanaka 
& Shibazaki 1996 and Campana et al. 1998 and references 
therein). Type I X-ray bursts were discovered during the de- 
clining phase of an outburst (Koyama et al. 1981), testifying 
to the presence of a neutron star. Peak X-ray luminosities are 
in the ~ (1 -4) x 10 37 ergs" 1 range (for the ~ 2.5 kpc distance 
inferred from its optical counterpart; Thorstensen et al. 1978). 
Close to the outburst maximum the X-ray spectrum is soft with 
an equivalent bremsstrahlung temperature of k 7t, r ~ 4 - 5 keV. 



Sporadic observations of Aql X-l during the early part of the 
outburst decay (Czerny et al. 1987; Tanaka & Shibazaki 1996; 
Verbunt et al. 1994) showed that when the source luminosity 
drops below ~ 10 36 ergs" 1 the spectrum changes to a power 
law with a photon index of T ~ 2, extending up to energies of 
- 100 keV (Harmon et al. 1996). ROSAT PSPC observations 
revealed Aql X-l in quiescence on three occasions at a level of 
~ 10 33 ergs" 1 (0.4-2.4 keV; Verbunt et al. 1994). In this lower 
energy band the spectrum is considerably softer and consistent 
with a black body temperature of kT^ ^0.3 keV. 

2. X-RAY OBSERVATIONS 

An outburst from Aql X-l reaching a peak luminosity of 
~ 10 37 ergs" 1 (2-10 keV) was discovered and monitored start- 
ing from mid-February, 1997 with the RossiXTE All Sky Mon- 
itor (ASM; Levine et al. 1997). Six observations were carried 
out with the BeppoSAX Narrow Field Instruments (NFIs) start- 
ing from March 8' /! , 1997 (seeTable 1), with the aim of studying 
the final stages of the outburst decay. Fig. la shows the light 
curve of the Aql X-l outburst as observed by the RossiXTE 
ASM and BeppoSAX MECS. The first part of the outburst can 
be fit by a Gaussian with sigma ~ 17 d. This is not uncom- 
mon in SXRTs (e.g. in the case of 4U 1608-52; Lochner & 
Roussel-Dupre 1994). 

The flux decay rate changed dramatically around MJD 505 12 
(March 5 ,/! , 1997). At the time of the first BeppoSAX observa- 
tion (which started on March 8''\ 1997) the source luminos- 
ity was decreasing very rapidly, fading by about 30% in 1 1 hr, 
from a maximum level of ~ 10 35 ergs" 1 . The second observa- 
tion took place on March 12' /! , 1997 when the source, a factor 
of ~ 50 fainter on average, reduced its flux by about 25% in 12 
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Fig. 1. — Light curve of the Feb. -Mar. 1997 outburst of Aql X-l (panel a. Data before and after MJD 50514 were collected with 
the RossiXTE ASM (2-10 keV) and the BeppoSAX MECS (1 .5-10 keV), respectively. RossiXTE ASM count rates are converted to 
(unabsorbed) luminosities using a conversion factor of 4 x 10 35 ergs" 1 (before MJD 50512) and 2 x 10 35 ergs" 1 (after MJD 50512) 
as derived from RossiXTE spectral fits (Zhang, Yu & Zhang 1998). BeppoSAX luminosities are derived directly from the spectral 
data (see text). The evolution of the flux from MJD 50480 to MJD 50512 is well fit by a Gaussian centered on MJD 50483.2. This fit 
however does not provide an acceptable description for later times (see the dot-dashed line), not even if the accretion luminosity is 
calculated in the propeller regime (dashed line). The straight solid line represents the X-ray luminosity corresponding to the closure 
of the centrifugal barrier L min (for a magnetic field of 10 8 G and a spin period of 1.8 ms) and the straight dashed line the luminosity 
gap due to the action of the centrifugal barrier, L cor . The dotted line marks the minimum luminosity in the propeller regime (Li c ). 
Panel b shows the BeppoSAX unfolded spectra of Aql X-l during the early stages of the fast decline (1) and during the quiescent 
phase (3-6, summed). The best fit spectral model (black body plus power law) is superposed to the data. 



hr. In the subsequent four observations the source luminosity 
attained a constant level of ~ 6 x 10 32 ergs" 1 , consistent with 
previous measurements of the quiescent luminosity of Aql X-l 
(Verbunt et al. 1994). The sharp decrease after MJD 50512 is 
well described by an exponential decay with an e-folding time 
~1.2d. 

The BeppoSAX LECS, MECS and PDS spectra during the 
fast decay phase, as well as those obtained by summing up 
all the observations pertaining to quiescence, can be fit with a 
model consisting of a black body plus a power law (see Table 2 
and Fig. lb). The soft black body component remained nearly 
constant in temperature (kTy^ ~ 0.3-0.4 keV), but its radius de- 
creased by a factor of ~ 3 from the decay phase to quiescence. 
The equivalent radius in quiescence (/?bb ~ 1 km) is consistent 
with the ROSAT results (Verbunt et al. 1994). The power law 
component changed substantially from the decay phase to qui- 
escence: during the decay the photon index was r <~ 2, while in 
quiescence it hardened to T ~ 1 . The two values are different 
with > 90% confidence (see Table 1). The ratio of the 0.5-10 
keV luminosities in the power law and black body components 
decreased by a factor of five between the first BeppoSAX ob- 
servation and quiescence. 



3. DISCUSSION 

The BeppoSAX observations enabled us to follow for the 
first time the evolution of a SXRT outburst down to quiescence. 
The sharp flux decay leading to the quiescent state of Aql X-l is 
reminiscent of the final evolution of dwarf novae outbursts (e.g. 
Ponman et al. 1995; Osaki 1996), although there are obvious 
differences with respect to the X-ray luminosities and spectra 
involved in the two cases, likely resulting from the different 
efficiencies in the gravitational energy release between white 
dwarfs and neutron stars. 



Models of low mass X-ray transient outbursts hosting an 
old neutron star or a black hole are largely built in analogy 
with dwarf novae outbursts. In particular, van Paradijs (1996) 
showed that the different range of time-averaged mass accretion 
rates over which the dwarf nova and low mass X-ray transient 
outbursts were observed to take place is well explained by the 
higher level of disk irradiation caused by the higher accretion 
efficiency of neutron stars and black holes. However, the out- 
burst evolution of low mass X-ray transients presents impor- 
tant differences. In particular, the steepening in the X-ray flux 
decrease of Aql X-l has no clear parallel in low mass X-ray 
transients containing Black Hole Candidates (BHCs). The best 
sampled light curves of these sources show an exponential-like 
decay (sometimes with a superposed secondary outburst) with 
an e-folding time of <~ 30 d and extending up to four decades 
in flux, with no indication of a sudden steepening (Chen et al. 
1997). In addition, BHC transients display a larger luminosity 
range between outburst peak and quiescence than neutron star 
SXRTs (Garcia et al. 1998 and references therein). Being the 
mass donor stars and the binary parameters quite similar in the 
two cases, it appears natural to attribute these differences to the 
different nature of the underlying object: neutron stars possess 
a surface and, likely, a magnetosphere, while BHCs do not. 

When in outburst accretion down to the neutron star surface 
takes place in SXRTs, as testified by the similarity of their prop- 
erties with those of persistent LMXRBs, especially the occur- 
rence of type I bursts and the X-ray spectra and luminosities. 
The mass inflow rate during the outburst decay decreases, caus- 
ing the expansion of the magnetospheric radius, r m . Thus, ac- 
cretion onto the neutron star surface can continue as long as 
the centrifugal drag exerted by the corotating magnetosphere 
on the accreting material is weaker than gravity (Illarionov & 
Sunyaev 1975; Stella et al. 1986). This occurs above a lim- 
iting luminosity L min = GMM min /R - 4 x 10 36 BjP^ 3 ergs" 1 , 
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Table 1 



Obs./Date 


LECS/MECS-PDS 


LECS Count Rate 


MECS Count Rate 


PDS Count Rate 




Expos, (s) 


(c s" 1 ) 


(c s" 1 ) 


(c s" 1 ) 


1/March 8"\ 1997 


5240/21342 


0.84 ±0.02 


2.2±0.01 


0.87 ±0.06 


2/March 12'*, 1997 


3247/21225 


(2.5 ± 0.4) x 10" 2 


(7.4 ±0.2) x 10~ 2 


< 0.19 


3/March 17'*, 1997 


5755/17258 


(5.6 ± 1.7) x 10" 3 


(1.3 ±0.1) x 10" 2 


< 0.24 


4/March20" 1 , 1997 


4287/22589 


(5.8 ± 1.9) x 10" 3 


(1.6±0.1) x 10" 2 


< 0.17 


5/April2" d , 1997 


8440/23576 


(6.2 ± 1.3) x 10" 3 


(1.2±0.1) x 10" 2 


< 0.17 


6°/May 6 th , 1997 


11789/21703 


(6.7±l.l)x 10" 3 




< 0.19 



° No MECS data were obtained. 



Table 2 

Summary of spectral fits. 

Obs.* Black body Black body Power law PL/BB Mean Luminosity t Reduced 
fcr bb (keV) flbb(km) T flux ratio* (erg s" 1 ) X 2 

1 0.42 ±0.02 2.6±0.3 1.9±0.1 3J 9 x 10 34 L0 

2 0.3^ 1.8±0.7 1.6 2 x 10 33 0.9 
3-6 0.3 ±0.1 0.8ffi 1.0 ±0.3 0.7 6 x 10 32 1.3 



Errors are la. 

* Spectra from the LECS and MECS (and PDS for the first observation) detectors have been considered. The spectra corresponding 
to the quiescent state have been summed up, in order to increase the statistics and an upper limit from the summed PDS data was 
also used to better constrain the power law slope. 

t Unabsorbed X-ray luminosity in the energy range 0.5-10 keV. In the case of the first observation the PDS data were included in 
the fit (the unabsorbed 0.5-100 keV luminosity amounts to 2 x 10 35 ergs" 1 ). 
\ Power law to black body flux ratio in the 0.5-10 keV energy range. 



where G is the gravitational constant; M, R, B = B$ 10 G and 
P = P-3 10" 3 ms are the neutron star mass, radius, magnetic field 
and spin period, respectively (here and in the following we as- 
sume M = 1.4 Mq and R = 10 6 cm). As r m reaches the coro- 
tation radius, r COI , accretion onto the surface is inhibited and a 
lower accretion luminosity (< L m j n ) of L cor = GMM m ; n / 

'"cor rs ~ J 

2 x 10 36 BgPl3 ergs" 1 is released. After this luminosity gap 
the source enters the propeller regime. If the mass inflow rate 
decreases further, the expansion of r m can continue up to the 
light cylinder radius, r\ c , providing a lower limit to the ac- 
cretion luminosity that can be emitted in the propeller regime 

Lie = GMMic Aic ~ 2 x 1 34 B\ P^ /2 erg s" 1 . Below L k the radio 
pulsar mechanism may turn on and the pulsar relativistic wind 
interacts with the incoming matter pushing it outwards. Matter 
inflowing through the Roche lobe is stopped by the radio pul- 
sar radiation pressure, giving rise to a shock front (Illarionov & 
Sunyaev 1975; Shaham & Tavani 1991). Clearly these regimes 
have no equivalent in the case of black hole accretion. 

3.1. The onset of the propeller 

During the February-March 1997 outburst of Aql X-l, 
RossiXTE observations led to the discovery of a nearly co- 
herent modulation at ~ 550 Hz (~ 1.8 ms) during a type I 
X-ray burst. A single QPO peak, with a centroid frequency 
ranging from uq P0 ~ 750 to 830 Hz, was also observed at 
two different flux levels, when the persistent luminosity was 
- 1.2 x 10 36 ergs" 1 and - 1.7 x 10 36 ergs" 1 (Zhang etal. 1998). 
In the presence of a single QPO peak, the magnetospheric and 
sonic point beat frequency models (Alpar & Shaham 1985; 
Miller et al. 1997) interpretation is ambiguous in that the QPO 



peak could represent either the Keplerian frequency at the in- 
ner disk boundary or the beat frequency. Moreover, the burst 
periodicity at ~ 550 Hz may represent the neutron star spin 
frequency, v s , or half its value (Zhang et al. 1998). In either 
cases, the possibility that accretion onto the neutron star surface 
takes place even in the quiescent state of Aql X- 1 faces serious 
difficulties: for a quiescent luminosity of order 10 33 ergs" 1 a 
magnetic field of only ^ 5 x 10 6 G would be required, in or- 
der to fulfill the condition r m £ r COI . For such a low magnetic 
field, Aql X-l and, by inference, LMXRBs with kHz QPOs can 
hardly be the progenitors of recycled MSPs. More crucially, the 
marked steepening in the outburst decay that takes place below 
~ 1 x 10 36 ergs" 1 , is accompanied by a marked spectral hard- 
ening, resulting from a sudden decrease of the flux in the black 
body spectral component. This is clearly suggestive of a tran- 
sition taking place deep in the gravitational well of the neutron 
star, where most of the X-rays are produced. The most appeal- 
ing mechanism is a transition to the propeller regime, where 
most of the inflowing matter is stopped at the magnetospheric 
boundary (Zhang, Yu & Zhang 1998). In Fig. la, the luminosity 
at MJD 50512 is identified with L min and the lower horizontal 
lines indicate the luminosity interval during which Aql X-l is 
likely in the propeller regime. 

Additional information on the neutron star magnetic field 
(and spin) can be inferred as follows. The observed ratio of 
the luminosity, Lqpo, when the QPO at ~ 800 Hz were de- 
tected and the luminosity L m ; n when the centrifugal barrier 
closes is LQ P0 /L m i n <~ 1.2-1.7. At L m i n the Keplerian fre- 
quency of matter at the magnetospheric boundary is, by defi- 
nition, equal to the spin frequency, i.e. v s ~ 550 or ~ 275 Hz 
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for Aql X-l. Based on beat-frequency models, at Lqpo the 
Keplerian frequency at the inner disk boundary can be ei- 
ther vk,qpo ~ 800 Hz or vk.qpo ~ v s + 800 Hz, depending on 
whether the single kHz QPOs observed corresponds to the Ke- 
plerian or the beat frequency. In the magnetospheric beat- 
frequency models, simple theory predicts that the Keplerian 
frequency at the magnetospheric boundary is oc L 3/ " 7 ; in the ra- 
diation pressure-dominated regime relevant to the case at hand, 
the Ghosh and Lamb (1992) model predicts instead oc L 3 / 13 . 
Therefore we expect vk,qpo/v s ~ {Lqpo/L^) 3 ! 1 ~ 1.2 and 

vk,qpoI v s ~ (LQpo/L m in) 3 / 13 ~ 1.1, in the two models, respec- 
tively. Such a low ratio clearly favors the interpretation in which 
vk.qpo ~ 800 Hz and v s <~ 550 Hz. In the sonic point beat- 
frequency model (Miller et al. 1997), the innermost disk radius 
is well within the magnetosphere, implying that the Keplerian 
frequency at the magnetospheric boundary is < vk.qpo- ln this 
case all possible combinations of v s and i>k,qpo are allowed. 
By using the observed L m ; n , a neutron star magnetic field of 
B ~ 1-3 x 10 8 G (depending on the adopted model of the disk- 
magnetosphere interaction) is obtained in the case v s ~ 550 Hz 
and B ~ 2-6 x 10 8 G in the case v s ~ 275 Hz. 

Once in the propeller regime, the accretion efficiency de- 
creases further as the magnetosphere expands for decreasing 
mass inflow rates (L acc oc M 9//? ). The ~ 1 d exponential-like 
luminosity decline observed with BeppoSAX is considerably 
faster than the propeller accretion luminosity extrapolated from 
the first part of the outburst (e.g. the Gaussian profile shown 
by the dashed line in Fig. la). We note here that the spec- 
tral transition accompanying the onset of the centrifugal barrier 
may also modify the irradiation properties of the accretion disk, 
contributing to X-ray luminosity turn off. Alternatively, an ac- 
tive contribution of the "propeller" mechanism or the neutron 
star spin-down energy dissipated into the inflowing matter can- 
not be excluded. 

3.2. A turned-on rotation-powered pulsar? 

It is unlikely that the quiescent luminosity of Aql X-l is pow- 
ered by magnetospheric accretion in the propeller regime. As 
shown in Fig. la, the quiescent X-ray luminosity is proba- 
bly lower than the minimum magnetospheric accretion lumi- 
nosity Lie allowed in the propeller phase (this remains true for 
B > 6 x 10 7 G if v s ~ 550 Hz, and for B > 3 x 10 8 G if v s ~ 275 
Hz). Moreover the BeppoSAX X-ray spectrum shows a pro- 
nounced decrease in the power law to black body flux ratio to- 
gether with a flattening of the power law component between 



the fast decay phase and quiescence, suggesting that a transition 
to shock emission from the interaction of a radio pulsar wind 
with the matter outflowing from the companion star has taken 
place. Note that an X-ray spectrum with a slope of T ~ 1.5 has 
been observed from the radio pulsar PSR B 1259-63 immersed 
in the wind of its Be star companion. Models of this interaction 
predict that a power law X-ray spectrum with a slope around 
r ~ 1.5 should be produced for a wide range of parameters 
(Tavani & Arons 1997). 

The additional soft X-ray component observed during the 
outburst decay (see Table 2) might be emitted at the polar caps 
as a result of the residual neutron star accretion in the pro- 
peller phase. Note that the equivalent black body radius de- 
creases for decreasing X-ray luminosities, just as it would be 
expected if the magnetospheric boundary expanded. Alterna- 
tively, the black body-like spectral component observed in qui- 
escence could be due to the streaming of energetic particles that 
hit the polar caps, in close analogy to the soft X-ray component 
observed, in MSPs, at the weaker level of ~ 10 30 - 10 31 ergs" 1 
(Becker & Triimper 1997). 

Assuming a magnetic field in the range derived in section 3.1 
(i.e. B ~ 1-3 x 10 8 G for u s = 550 Hz and B ~ 2-6 x 10 8 G 
for v s = 275 Hz), we can consistently explain the <~ 10 33 ergs -1 
quiescent X-ray luminosity as powered by a radio pulsar en- 
shrouded by matter outflowing from the companion star, if 
the conversion efficiency of spin-down luminosity to X-ray is 
~ 0. 1 - 10%. This is consistent with modeling and observations 
of enshrouded pulsars (Tavani 1991; Verbunt et al. 1996). 

There are chances of observing a MSP (a simple scaling from 
MSPs implies a signal at 400 MHz of ~ 10 mJy; see Kulkarni 
et al. 1990), even though the emission would probably be spo- 
radic, like in the case of the pulsar PSR 1744-24A due to the 
large amount of circumstellar matter (see Lyne et al. 1 99 1 ; Sha- 
ham& Tavani 1991). 

In summary Aql X-l appears to provide the first example of 
an old fast rotating neutron stars undergoing a transition to the 
propeller regime at first, followed by a transition to the radio 
pulsar regime, as the transient X-ray emission approaches its 
quiescent level. Therefore, Aql X-l (and possibly SXRTs in 
general) likely represents the long-sought "missing link" be- 
tween LMXRBs and recycled MSPs. 

We thank J. van Paradijs, FK. Lamb and an anonymous ref- 
eree for useful comments. S.N. Zhang, W. Yu and W. Zhang 
shared with us their RossiXTE results in advance of publica- 
tion. 
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